Magnetic Excitations in the High T c Iron Pnictides 
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We calculate the expected finite frequency neutron scattering intensity based on the two-sublattice 
collinear antiferromagnet found by recent neutron scattering experiments as well as by theoretical 
analysis on the iron oxypnictide LaOFeAs. We consider two types of superexchange couplings 
between Fe atoms: nearest-neighbor coupling J\ and next-nearest-neighbor coupling Ji. We show 
how to distinguish experimentally between ferromagnetic and antiferromagnetic J\. Whereas mag- 
netic excitations in the cuprates display a so-called resonance peak at (n, n) (corresponding to a 
saddlepoint in the magnetic spectrum) which is at a wavevector that is at least close to nesting 
Fermi-surface-like structures, no such corresponding excitations exist in the iron pnictides. Rather, 
we find saddlepoints near (ir,n/2) and (0,7r/2) (and symmetry related points), which are not close 
to nesting the Fermi surfaces. 

PACS numbers: 74.25.Ha, 74.70.-b, 75.30.Ds, 76.50.+g 
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The recent discovery of superconductivity exceeding 
50K in a new class of materials holds tremendous po- 
tential for understandingthe origin of high temperature 
superconductivity, [l], [3,11, 0, H, Q Similar to the cuprate 
superconductors, the iron pnictides also have a layered 
structure, and display magnetism in the undoped parent 
compound. Both become superconducting upon doping. 
And like the cuprates, the transition metal layer is be- 
lieved to play an important role in the superconducting 
pairing. On the other hand, the parent compound of 
LaOFeAs is a poor metal at room temperature, as op- 
posed to a correlated insulator as in the cuprates. 

Initially band structure calculations suggested the ma- 
terials are nonmagnetic but close to a strong magnetic 
instability. @, Bill However, subsequent calculations have 
shown that the antiferromagnetic state has lower energy 
than the nonmagnetic state because of Fermi surface 
nesting. hJ, U, HJ In Ref. 12, a stripe-like antiferromag- 
netic ground state was suggested based on strong nesting 
effects. Recent neutron scattering experiments [13j have 
shown that the parent compound of LaOFeAs is a long- 
range ordered antiferromagnet with a type of spin stripe 
order (i.e. unidirectional spin density wave). However 
the magnetic moment was found to be 0.36(5)/ib per 
iron, which is much smaller than the calculated value of 
- 2.3^ B per iron. [13, El, E3] ■ 



From an analysis of the superexchange interactions, 
Ref. llJ suggested that the next-nearest-neighbor in- 
teraction J2 is antiferromagnetic (AFM), while the 
nearest- neighbor interaction J\ is ferromagnetic (FM). 
However a first-principles band structure calculation 
predicts that the nearest-neighbor interaction is also 
antiferromagnetic. [lil flij ] They predict that | J2I is al- 
most as twice large as J\. In both cases, the competition 
between J\ and J2 leads to a type of stripe-ordered two- 
sublattice antiferromagnetic ground state (Fig P) when 
I J2/J1I is larger than the critical value. 17|,[l8( While the 
interactions J\ and J2 can compete, the uniaxial SDW 




FIG. 1: (Color online) Two-sublattice collinear antiferromag- 
net on the Fe-square lattice. Shaded region is the magnetic 
unit cell. 



considered in Fig. Q] is a classical ground state of the sys- 
tem, and it is thus not frustrated in the sense of having 
a macroscopic ground state degeneracy. 

We use linearized spin wave theory to calculate the 
magnetic excitations and sublattice magnetization for the 
two-sublattice collinear antiferromagnet with nearest- 
neighbor superexchange coupling J\ and antiferromag- 
netic next-nearest- neighbor superexchange coupling J2. 
We present results for ferromagnetic coupling J\ as well 
as for antiferromagnetic coupling J\. (See Fig. [TJ) We 
find the results are quite different for the two cases, so 
that comparing our calculations with future neutron scat- 
tering results at finite frequency will be able to distin- 
guish these two cases. 

The model Hamiltonian is described by the Heisenberg 
spin model on the square lattice 
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where < ij > nn and < ij > nnn mean the nearest- 
neighbor and next-nearest-neighbor spin pairs respec- 
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tively. There are two spins in each unit cell, as shown in 
Fig. [TJ We study the elementary excitations of the clas- 
sical ground state of this model by using the well-known 
Holstein- Primakoff boson method. The dispersion and 
intensities are calculated by quantizing the classical spin 
waves. 

We use Holstcin-Primakoff bosons to quantize about 
the collinear antiferromagnetic ground state found in re- 
cent neutron scattering. [Ij] 
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H = Eci + S^JA^a+a^ + -(S k a+a+ k + Bl k a k a_ k ] 

k 

(2) 

where Eci = —2J2NS 2 is the classical ground state en- 
ergy and 

A* = (4J 2 +2JicosA x ), (3) 
£? k = (2 Ji cos k y + 4J 2 cos k x cos k y ). (4) 

We can diagonalize the Hamiltonian using the Bogoli- 
ubov transformation 



fe k = cosh# k a k — sinh0 k a_ k . 
The diagonalized Hamiltonian is 



H 



&k + Eci + Eq 



where w(k) is the spin wave dispersion 



w(k) - SJAi - B 



and Eo is the quantum zero-point energy correction 

£o = f]T(-A k + c(k)). 

k 

For |Ji| = 1 and J 2 = 2, we get E = -0.3327V5. 
We find that there is only one spin wave band 

u(k x , k y ) = 



(5) 



(6) 
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(8) 



2Sy (2 J2 + Ji cos k x ) 2 — ( J\ cos k y + 2 J2 cos k x cos k y ) 2 
The associated spin wave velocities are 



2SVi + 2J 2 



4/2 



(10) 

(11) 



Notice that v x becomes imaginary for | J\\ > 2| J%\, indi- 
cating a change in the classical ground state configura- 
tion. 

Fig.[2]shows the spin wave band with the nearest neigh- 



bor coupling both antiferromagnetic (Fig. 2(a)) and fer- 
romagnetic (Fig. 2(b) ). The presence of saddlepoints can 
be seen, and we will return to this point later. In addi- 
tion, because the (-7T, it) point is a magnetic reciprocal 
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FIG. 2: (Color online) Spin-wave dispersion band for the two- 
sublattice collinear antiferromagnet shown in Fig. [1] (a) Dis- 
persion with both couplings antiferromagnetic. Here we have 
set Ji = 1 (AFM) with J 2 = 2 (AFM). (b) Dispersion with 
ferromagnetic nearest neighbor coupling. Here we have set 
Ji = -1 (FM) with ,h = 2 (AFM). 



lattice vector, the dispersion must have uj — > at this 
point, although as we will see there is no zero- frequency 
(9) intensity associated with this part of the dispersion. This 
precludes finite frequency weight at the (tt, it) point from 
this band. 

We calculate the zero-temperature dynamic structure 
factor using the same method, 0,11 



S(k,u;) = £ £ \<f\S\k)\Q>\ 2 6(u-Lo f ). (12) 

/ i=x,y,z 

Here |0 > is the magnon vacuum state and |/ > denotes 
the final state of the spin system with excitation energy 
u)f. S z does not change the number of magnons, con- 
tributing to the elastic part of the structure factor. S x (k) 
and S y (k) contribute to the inelastic dynamic structure 
factor through single magnon excitations. 

In Figs. [3] and 21 we show the expected neutron scat- 
tering intensity for constant energy cuts in k-space. We 
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FIG. 3: (Color online) Constant-energy cuts (untwinned) of 
the dynamic structure factor S(k,ui) for J2 = 2 (AFM): (a) 
Ji = -1 (FM), (b) Ji = 1 (AFM). The x-axis and y-axis 
correspond to k x and k y respectively with the range (0, 2-n). 
We have integrated over an energy window of ±0.2\Ji\S. 



FIG. 4: (Color online) Constant-energy cuts (twinned) of 
the dynamic structure factor S(k, u>) for .1-2 = 1 (AFM): (a) 
Ji = -1 (FM), (b) Ji = 1 (AFM). The x-axis and y-axis 
correspond to k x and k y respectively with the range (0, 27r). 
We have integrated over an energy window of ±0.2|Ji|S. 



show our predictions from spin wave theory for both fer- 
romagnetic and antiferromagnetic J\. Fig. [3] shows the 
expected neutron scattering intensity from a single do- 
main of the magnetic order {i.e. for an untwinned case), 
and Fig. |4] shows the expected scattering intensity for the 
case where there is an equal contribution from domains 
with both orientations of the magnetic order (i.e. for a 
twinned case). 

For ferromagnetic J\, at low frequency, the strongest 
diffraction peaks are located at (0, tt). (See Fig. [3]) How- 
ever more intensity weight shifts to (tt, 0) when J\ is an- 
tiferromagnetic. There is also a spin wave cone emerging 
from (tt, 7r), but the intensity is much weaker than the 
cones emanating from other magnetic reciprocal lattice 
vectors, since zero frequency weight is forbidden at (tt, tt) 
for the magnetic order we consider. At high energy, the 
difference between ferromagnetic J\ and antiferromag- 
netic Ji becomes more apparent. For example, for FM 
Ji, there are two strong spots along the (tt, k y ) direction, 
whereas for AFM J%, they are along the (0, k y ) direction. 
In real materials, stripe order can be twinned due to, 
e.g., a finite correlation length, local disorder pinning, or 
crystal twinning. Therefore we show the twinned con- 
stant energy cut plots in Fig. [4] for both FM and AFM 
nearest neighbor coupling J\ . 

As can be seen from the dispersion in Fig. [2j there 
are saddlepoints in the spin wave excitation spectrum at 
various points in fc-space. For the case of both couplings 
antiferromagnetic, these occur at (tt/2,0) and (tt,tt/2) 
and symmetry related points. For ferromagnetic nearest 
neighbor coupling, saddlepoints can be seen at (0,n/2) 
along with weak saddlepoints possible at (tt/2,0) and 
(7r/2,7r) and symmetry related points. The integrated 
intensity is generally large at such saddlepoints. In the 



cuprates, there is a saddlepoint in the magnetic excita- 
tions at (tt, tt) which has been empirically connected to 
superconductivity, in that it increases in intensity at the 
onset of superconductivity, i. e. the "resonance peak" . 
There has been much discussion concerning this scatter- 
ing phenomenon in the cuprates, particularly because it 
is close to nesting vectors for the corresponding Fermi 
surface. However, in the case of the iron pnictides, the 
saddlepoints we find here are quite far from any nesting 
vectors. 

Experimentally, the magnetic moment per iron was 
found to be 0.36(5)/kb, which is much smaller than the 
expected value of ~ 2.3/is per iron site. 1(| U, 13] The 
zero point energy of the spin waves reduces the sublat- 
tice magnetization. It was suggested in Ref. that the 
competition between J\ and J2 may be responsible for 
the small moment observed in experiment. The sublat- 
tice magnetization m is defined as 



= < Sf >=S-Ar 



(13) 



where Am is the deviation of sublattice magnetization 
from the saturation value, 

Am = < af ai > 

= < a i ak > 

k 



1 v SA k 1 ^ SA k 1 

2V k ^ [ uj(k) Vk Y ^( k ) e/Mk) - 1 



(14) 



The first term f\ m i uantum comes from quantum zero 
point fluctuations. The sceond term Am ttermo ' comes 
from the classical thermal fluctuation, which is divergent 
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FIG. 5: IJ2/J1I dependence of the reduction of the sublattice 
magnetization due to zero point energy of the spin waves. 
Here we have used I Ji I = 1 



at any finite temperature in agreement with the Mermin- 
Wagner theorem. (The very presence of the broken sym- 
metry observed in experiment implies that there is some 
finite coupling between planes, however weak.) 



Here we calculate Am^ uantum by 



^ m quantum 



~2w~2n Lu(k) ~ 2' 



(15) 



It is difficult to get the analytical form of the integral. 
Thus we numerically calculate ^ m i uantum _ From the 
symmetry, the above integral does not change when Ji 



changes sign. In Fig. El Ar 



.quantum 



is plotted as a func- 



tion of the superexchange coupling ratio IJ2/J1I. It is 
S -independent. If S is in between 1 and |, it will re- 
duce the m by 13% - 20%. Am« nantMm decreases with 
increasing J2/J1 because stonger J 2 stablizes the two- 
sublattice collinear antiferromagnet state. This devia- 
tion is not sufficient to explain the observed value of the 
sublattice magnetization. 

In conclusion, we have used spin wave theory to cal- 
culate the magnetic excitations and sublattice magneti- 
zation for the two-sublattice collinear antifcrromagnetic 
state of the new La(0\- x F x )FeAs high-T c superconduc- 
tors. We have studied both ferromagnetic and antifer- 
romagnetic nearest-neighbor coupling J\ with antiferro- 
magnetic next- nearest- neighbor coupling Ji- We calcu- 
late the predicted inelastic neutron scattering pattern 
based on spin wave theory. Comparison with future in- 
elastic neutron scattering studies can be used to distin- 
guish the sign of J\. We find that the sublattice mag- 
netization can be reduced by the zero-point motion of 
spin waves, although not enough to account for the small 
moments observed in experiment. In addition, we iden- 



tify several saddlepoints in the magnetic excitation spec- 
trum. While magnetic excitations in these regions are 
expected to have extra intensity due to the saddlepoint 
structure, these corresponding wavevectors are not near 
nesting vectors of the Fermi surface. 
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Note added: Some results from spin wave calculations 
have also been reported by Ref. [2l| . 
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